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Using exact diagonalization (ED) and linear spin wave theory (LSWT) we study the influence of 
frustration and quantum fluctuations on the magnetic ordering in the ground state of the spin-i 
J1-J2 Heisenberg antiferromagnet {J1-J2 model) on the body-centered cubic (bcc) lattice. Contrary 
to the J1-J2 model on the square lattice, we flnd for the bcc lattice that frustration and quantum 
fluctuations do not lead to a quantum disordered phase for strong frustration. The results of both 
approaches (ED, LSWT) suggest a first order transition at J2/J1 ~ 0.7 from the two-sublattice Neel 
phase at low J2 to a coUinear phase at large J2. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The properties of the two-dimensional (2d) J1-J2 
model have attracted a great deal of interest during the 
last decade (see e.g. Refs. [lllllillllli,lllili[lfl,II3 

and references therein). The Hamiltonian of the J1-J2 
model is 

i/ = Ji^S,.S,+ J2^S,.Sj, (1) 

where Ji = 1 is the nearest-neighbor, and J2 > is the 
frustrating next-nearest-neighbor Heisenberg exchange. 
We are interested in the extreme quantum case, i.e. we 
consider the spin quantum number is s — 1/2. For the 
square lattice it seems to be well-accepted that there is 
a quantum spin-liquid phase between J2/J1 ~ 0.38 and 
J2/J1 « 0.60 and that the corresponding quantum phase 
transitions from the Neel ordered state to the spin-liquid 
state at J2/J1 ~ 0.38 is of second order. The nature of 
the transition from the spin-liquid state to the coUinear 
state at J2/J1 « 0.60 is still under discussion but there 
are indications that it might be of first order The 
Neel phase for small J2 is characterized by an antiparallel 
alignment of nearest-neighbor spins with a correspond- 
ing magnetic wave vector QncsI = (tt, tt). The coUinear 
state for large J2 is twofold degenerated and the corre- 
sponding magnetic wave vectors are Qj^j — (tt, 0) and 
Qco/ — (0,7r). The two coUinear states are character- 
ized by a parallel spin orientation of nearest neighbors 
in vertical (horizontal) direction and an antiparallel spin 
orientation of nearest neighbors in horizontal (vertical) 
direction and exhibit therefore Neel order within the ini- 
tial sublattices A and B. The properties of the spin- liquid 
phase are a current field of active research. Even addi- 
tional quantum phase transitions at (J2/J1) ~ 0.34 and 
(J2/J1) « 0.50 are discussed [Toj. 
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The properties of quantum spin systems strongly de- 
pend on the dimensionality. So, contrary to the 2d model, 
the one-dimensional Ji — J2 model does not have a Neel 
ordered ground state, but exhibits a transition from a 
critical state to a dimer phase at J2/J1 = 0.241 (see e.g. 
Clllllll)- Though the tendency to order IS more pro- 
nounced in three-dimensional (3d) quantum spin systems 
than in low-dimensional ones a spin-liquid phase is also 
observed for frustrated 3d systems like the Heisenberg 
antiferromagnet on the pyrochlore lattice 0, 0| . 

In this paper we consider the 3d version of the Ji — J2 
model. To cover the possibility to have in a 3d model 
both Neel phases found for the 2d model we need a 3d bi- 
bipartite lattice, i.e. a lattice consisting of two interpen- 
etrating bipartite sublattices. The bcc lattice consists of 
two interpenetrating, identical simple cubic sublattices. 
Each simple cubic sublattice consists of two interpene- 
trating, identical fee lattices. Therefore the bcc lattice is 
a 3d bi-bipartite cubic lattice. 

The classical ground state for the bcc lattice corre- 
sponds to that of the square lattice: For J2/J1 < (Xc it is 
a usual two-sublattice Neel state whereas for J2/ Ji > (Xc 
an antiferromagnet with four sublattices is realized. The 
transition point olc depends on the coordination num- 
bers etc = zi/(2z2), where z\ is the number of nearest- 
neighbors (Ji bonds) and Z2 of next-nearest-neighbors 
(J2 bonds). Consequently we have ac = 1/2 for the 
square lattice but etc = 2/3 for the bcc lattice. Therefore 
we define as the appropriate parameter of frustration 



and we have ^^9"'^'^'= = p^^^^ = 1/2. 

In what follows we use the exact diagonalization 
scheme and the linear spin wave theory to calculate the 
ground-state properties of the J\ — J2 model on the bcc 
lattice and compare the results with the corresponding 
ones for the square lattice. We will present the ground- 
state energy, the violation of the Marshall-Peierls sign 
rule, the sublattice magnetizations and the spin gap of 
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finite lattices in section II. Properties of the infinite lat- 
tices will be given in section III, where the results of the 
extrapolation of the ground-state energy, the sublattice 
magnetizations (from ED data) and the corresponding 
results of the LSWT are shown. 



II. EXACT DIAGONALIZATION 

A. The generation of finite bcc lattices 

The generation of finite 3d lattices with periodic 
boundary conditions is less transparent than for 2d lat- 
tices. As has been recently pointed out by Betts and 
coworkers 0, 0, 0| the use of a triple of edge vec- 
tors in upper triangular lattice form (utlf) leads to 
a systematic generation of finite 3d lattices. In this pa- 
per we use the utlf edge vectors and follow strictly Refs. 
[isl Il9l | . Finite parallelepipeds that build up the infinite 
bcc lattice can be defined by three edge vectors, 

3 

La = ^ riafi^p, (3) 

13=1 

where with (3 = 1,2,3 are integers and ai = 

(1,1,-1), sl2 = (1,-1,1), a3 = (—1,1,1) are the ba- 
sis vectors of the lattice connecting nearest neighbors. 
The lattice vectors connecting next nearest neighbors are 
bi = (±2,0,0),b2 = (±0,2,0),b3 = (0,0, ±2). 

There are altogether 10 finite bcc lattices with iV < 36 
listed in Tabled which fulfill the following three condi- 
tions: (i) Every site i of the bcc lattices should have 8 
nearest and 6 next-nearest neighbors, which means that 
they have the full number of nearest and next-nearest 
neighbors, (ii) The finite lattices should be bi-bipartit 
in order to avoid frustration due to boundary conditions 
for p = and p — > oo. (Notice, that finite bcc lattices 
may be not bi-bipartite even though the infinite bcc lat- 
tice is.) (iii) Furthermore they should be topologically 
distinct, i.e. the spin Hamiltonian should exhibit dif- 
ferent physical properties. 

B. Ground-state energy 

The ground-state energy gives first insight in the na- 
ture of possible zero-temperature phase transitions. In 
the thermodynamic limit a kink in Eq{p) (respectively a 
jump in the first derivative dEo/dp) signals a first order 
transition, whereas a smooth dEo/dp is compatible with 
second order transitions. Furthermore the maximum in 
Eq{p) indicates the point of maximal frustration. The 
ground-state energy for the classical model consists of 
two straight lines Eq°'^{p) — {p ~ 1)N for p < 0.5 and 
Eq°'^[p) — —pN for p > 0.5 with a kink (maximum) at 
p = 0.5. 

The quantum ground state |^'o) is a singlet eigen state 
of total spin for all J2 ■ In analogy to the square lattice 



TABLE I: The ten finite bcc lattices are used for exact diag- 
onafization. Li, L2, L3 are the three edge vectors in up- 
per triangle lattice form. is the Number of sites and 
7 — A, B ,C, D, ... is an additional label corresponding to a 
notation used in to distinguish finite lattices with identi- 
cal N. 
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N7 


Li 


L2 


L3 


24C 


(2,0,10) 


(0,2, 6) 


(0,0,24) 


28D 


(2,0,10) 


(0,2, 6) 


(0,0,28) 


32D 


(2,2, 4) 


(0,8, 0) 


(0,0, 8) 


32F 


(2,0, 6) 


(0,4, 8) 


(0,0,16) 


32H 


(2,0,10) 


(0,2, 6) 


(0,0,32) 


32J 


(4,0, 4) 


(0,4, 4) 


(0,0, 8) 


32K 


(2,0, 6) 


(0,4, 4) 


(0,0,16) 


36A 


(2,0,10) 


(0,2, 6) 


(0,0,36) 


36B 


(2,0,14) 


(0,2,10) 


(0,0,36) 


36C 


(2,2, 4) 


(0,6, 6) 


(0,0,12) 


TABLE II: 


The ground-state 


energy per site of the 10 con- 


sidered finite bcc lattices for different values of J2 (Ji = 1). 






Eo/N 




Nq 


J2 = 


J2 = 0.7 


J2 = 1.3333 


24C 


-1.21305 


-0.73883 


-1.34537 


28D 


-1.20223 


-0.73003 


-1.32544 


32D 


-1.19572 


-0.72851 


-1.31225 


32F 


-1.19512 


-0.72744 


-1.31217 


32H 


-1.19474 


-0.72708 


-1.31189 


32J 


-1.19440 


-0.74003 


-1.31688 


32K 


-1.19408 


-0.73203 


-1.31423 


36A 


-1.18953 


-0.72119 


-1.30021 


36B 


-1.19264 


-0.72180 


-1.30145 


36C 


-1.19278 


-0.72248 


-1.30149 



l^'o) of finite bcc lattices with N mod 8 = has 
the same translational symmetry for small and large J2 
(kgs = (0,0,0)), whereas the translational symmetry of 
I^E'o) for lattices with N mod 4 = but N mod 8 7^ is 
different for small and large J2 . The change of symmetry 
from kgs = (0,0,0) to k^s = (7r,7r,7r) appears slightly 
right from the maximum of the ground-state energy. 

For the quantum model we show exact-diagonalization 
results of Eq/N for three different values of frustration 
in Table ITU While J2 = and J2 = 1.3333 correspond 
to zero or small frustration, J2 = 0.7 is in the region of 
strong frustration. 

For the sake of clearness we present in the figures 
only results of selected finite lattices. To illustrate the 
finite-size effects in the most of the subsequent figures 
we present data for the smallest (N=24) and the largest 
lattices we have calculated, where for N = 32, 36 we have 
chosen the lattices having highest symmetry. Note, that 
the curves for lattices of identical N look very similar. 

For comparison with the square lattice we have recal- 
culated data of 0] up to iV = 6 X 6 = 36. However, we 
think a square lattice with = 4 x 4 = 16 is compa- 
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FIG. 1: Ground-state energy per site in units of Ji for the 
bcc lattices 24C, 36C and for the square lattices (scaled) with 
A'' = 16 (16sq) and A'^ = 36 (36sq) versus p. The inset is an 
enlargement of the strongly frustrated region around p = 0.52. 

rable to a bcc lattice with = 36. This can be seen if 
one looks at the characteristic lengths L^d oc N'^/^ and 
oc N-^/^. Further, we mention that the square lattice 
of TV = 16 contains 5 neighborhood shells, whereas some 
of the more dense finite bcc lattices with = 32 and 
TV = 36 contain even more. 

In Fig. Hone finds the ground-state energies of the bcc 
lattices 24C, 36C and for comparison the iV = 16 and the 
TV = 36 square lattices. To have comparable curves we 
scaled the ground-state energy of the square lattice with 
the factor 7/4. Fig.^ illustrates that, contrary to the 
2d model, the ground-state energy of the 3d bcc lattice 
behaves very similarly to the classical model. As can 
be seen in Fig.^ the kink in the ground-state energy is 
almost independent of the size of the bcc lattices. This 
can be interpreted as an indication that the kink survives 
in the thermodynamic limit. 

C. Ground-State Phase Relationships 

The phase relationships of the Ising basis states \n) in 
the ground state |^o) of the bipartite Heisenberg antifer- 
romagnet (i.e. J2 = in (J^)) follow the Marshall-Peierls 
sign rule "2ll| . This sign rule can be formulated as 

1*0) = Xl^"!'^)' c„>0. (4) 

n 

Here the Ising states \n) are defined by 

\n) = (-l)^/2-^'^(^)|mi} ® |to2) ® ... ® |mjv), (5) 

where jwi), i = l,...,iV, are the eigenstates of the site 
spin operator SI (i.e. = ±i) and M(X) = X^ifX "^i- 
The standard Marshall-Peierls sign rule appropriate for 
the Neel phase at small p is obtained for X — A, i.e. X 
labels one of the two equivalent sublattices. For large 



values of p we have antiferromagnetic order within the 
initial sublattices A and B and A and B resolve into 4 
sublattices {A — > A\, A2 and B — > Bi, B2). Then a 
modified sign rule holds with X — AiiJ Bi. 

As pointed out in 0, |^ and very recently in the 
sign rule may survive some frustration but is clearly vi- 
olated for the square lattice in the strongly frustrated 
spin-liquid region. Hence we can use the violation of the 
Marshall-Peierls sign rule as an indication of the break- 
down of the two-sublattice Neel state. 

In Fig.[21one finds the weight g{X) — Y^'„{cnY of the 
Ising states \n) fulfilling the Marshall-Peierls sign rule 
(i.e. the sum is restricted to the subset of states 
having c„ > 0) for two finite bcc lattices and two square 
lattices. For the Ji — J2 model on the bcc lattice the 
rule 101) is violated almost discontinuously at that point 
where the ground-state energy has its maximum. This 
is a hint to a very drastic change of the ground state on 
the bcc lattice around p = 0.52 which can be attributed 
to an abrupt breakdown of the two-sublattice Neel state. 
On the other hand for the Ji — J2 model on the square 
lattice the violation of the sign rule starts smoothly 
and becomes significant near p = 0.4 where the second- 
order transition to the spin-liquid state takes place. The 
modified sign rule with g{Ai U Bi) changes also discon- 
tinuously for the bcc lattice but smoothly for the square 
lattice. However, there is a violation of the modified rule 
also in the coUinear phase, which can be attributed to the 
coupling between the both antiferromagnetic subsystems 
living on the initial sublattices A and B. Only for large 
p ^ 1, where the ground state becomes a product state 
of both antiferromagnetic subsystems the modified rule 
is rigorously fulfilled. 




0.5 I ' ' ' 
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FIG. 2; The weight g{A) (that is 1 for p = 0) and g{Ax U Bi) 
(that is 1 for p ^ 1) of basis states |n) fulfilling the Marshall- 
Peierls sign rule for the bcc lattice 24C, 36C and the square 
lattice with A'^ = 16 (16sq) and A^ = 32 (32sq) versus p. 
On the bcc lattices the sign rule is completely fulfilled up to 
p = 0.4 and 99.8% fulfilled at p = 0.51. 
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TABLE III; Sublattice magnetizations and of the 10 
considered bcc lattices for different values of J2. 







mi 


Na 


J2 =0 


J2 = 1.3333 


24C 


0.09362 


0.11006 


28D 


0.09057 


0.10429 




U.UO ( ( 


u.uyyoo 


32F 


0.08800 


0.09959 


32H 


0.08811 


0.09964 


32J 


0.08819 


0.09897 


32K 


0.08827 


0.09933 


36A 


0.08600 


0.09605 


36B 


0.08516 


0.09562 


36C 


0.08509 


0.09561 



0.12 



0.08 















; 24C 

y ; 36C 

: V 16sq 

36sq — ^- - - 

^ 













0.2 0.4 0.6 0.8 
P 



FIG. 3: Sublattice magnetizations (maximal for small p) 
and (maximal for large p) of the bcc lattices 24C, 36C and 
the square lattices with A'^ = 16 (16sq) and TV = 36 (36sq) 
versus p. 
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FIG. 4: Sublattice magnetizations (maximal for small p) 
and jtIq (maximal for large p) of the bcc lattice 36C, 32D, 
28D and 24C near the phase transition point p ~ 0.52. 



D. Sublattice magnetizations 

Of course the most important parameter to study Neel 
ordering is the sublattice magnetization. In finite sys- 
tems the conventional antiferromagnetic LRO, that is re- 
alized for small p, has to be described by the square of 
the sublattice magnetization of one spin component 

with QATeei = (tt, tTjTt) for the bcc lattice (and Q^veei = 
(tt, tt) for the square lattice). For large values of p the 
magnetic wave vectors Q^l^^ — (±7r/2, ±7r/2, ±7r/2) have 
to be used for the bcc lattice (and Q;':^; = (tt, 0) or 
Qco/ = (O)'"') foi' the square lattice) to describe the 
coUinear phase with antiferromagnetic order within the 
initial sublattices A and B. We denote in the following 
the order parameter of the Neel phase calculated with 
Qwee; with and that of the coUinear phase calculated 
with Q;!:^; or Q^^j with m^. Notice, that is identi- 
cal for Q;':^; and Q^^; . In Table IIIII we give the order 
parameters of the ten finite bcc lattices for J2 = and 
J2 = 1.3333. The behavior of rv? and shown in Fig.|2| 
again illustrates, that the influence of the frustration on 
the ground-state properties is basically different for the 
square and bcc lattice and suggests a direct first-order 
transition between both Neel phases for the bcc lattice. 
A more detailed presentation of the transition region is 
given in Fig.^] One finds that the position of the tran- 
sition only slightly depends on size and symmetry of the 
finite lattices. Moreover, the width of transition region is 
getting smaller with growing N . Again we mention that 
the region of transition is related to the maximum of Eq 
and the significant violation of the Marshall-Peierls sign 
rule. 
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FIG. 5: Spin gap Ast in units of Ji, i.e. the gap to the first 
triplet excitation for the bcc lattices 24C, 36C and the square 
lattices with N=16 (16sq) and N=32 (32sq) versus p. 
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E. Spin gap 

Another indication for a possible quantum disordered 
spin-liquid state is the spin gap, i.e. the gap Ast between 
the singlet ground state and the first triplet excitation. 
In Neel ordered systems we have Goldstone modes and no 
spin gap is observed in the thermodynamic limit. Con- 
trary to this, quantum disorder is accompanied by the 
opening of a spin gap. We show the gap Ast of two 
finite bcc lattices in Fig.|31 The first triplet excitation 
relevant for the gap belongs to the translational symme- 
try kt = (tt, tt, it) for small p and kt = (ti'/S, 7r/2, tt/2) for 
large p. For comparison we show the gap for the square 
lattice of -/V = 16 and = 32 sites. Of course, the gap 
of a finite lattice is finite. However, in the long-range 
ordered Neel and coUinear phases the extrapolation to 
the thermodynamic limit yields a vanishing gap. Obvi- 
ously there is no increase in the gap for the bcc lattices 
till about p « 0.52 where the transition from the two- 
sublattice Neel phase to the coUinear phase takes place. 
In the coUinear phase the relevant coupling parameter 
for excitations is J2 instead of Ji and consequently the 
gap increases linearly with p. Clearly we see that there 
is no special behavior of Ast near the transition point 
p « 0.52. The gap for the square lattice shows a simi- 
lar behavior for parameter regions where magnetic long- 
range order in the ground state is present, but around 
p = 0.5 the behavior of Ast is in contrast to the bcc lat- 
tice. In this region a quantum disordered gapped phase 
for the square lattice is expected which is consistent with 
the significant increase of Ast near p = 0.5. 

We conclude from the examination of the spin gap that 
there are no indications for a quantum disordered gapped 
phase for the bcc lattice. 



III. INFINITE BCC LATTICES 



The same relation is valid for m^. The results are 
presented in Figs .[HI and [7| The discussion of the data is 
given below. 

B. Linear spin wave theory (LSWT) 

1. LSWT for small J2 

Starting from the classical two-sublattice Neel state 
we choose a two-boson representation of the Hamilto- 
nian (^). Rewriting the spin operators of the Hamil- 
tonian in terms of bose operators by using the usual 
Holstein-Primakoff transformation and taking into ac- 
count only quadratic terms in the bose operators, we 
obtain a bosonic Hamiltonian in Fourier transformed rep- 
resentation 

H = s'A^(-8 + 6p) + ^{Ak(a^ak + 6^6k) (9) 

k 

-HBk(ak&-k + <^'ik)} 
with the coefficients 

= .5(8-6p(l-72k)) (10) 
Bk = 8571k. (11) 

The structure factors of nearest neighbors and next- 
nearest neighbors are given by 

7ik — cos kx cos ky cos kz (12) 

72k = -(cos2fca; -f cos2A;j^ + cos2fc2). (13) 
o 

The ground-state energy per site is then 

So/A^ = s'(-8 + 6p) + ^^(c^k-^k) (14) 

k 



A. Finite-Size Extrapolation 

To obtain properties of the infinite bcc lattice we ex- 
trapolate the ED data of all ten lattices listed in Tabled 
The finite-size extrapolation is a well elaborated approx- 
imation scheme successfully applied to many 2d quan- 
tum spin systems like the Ji — J2 model on the square 
lattice ,2j. But even for 3d lattices this scheme may 
lead to precise data for the infinite lattice [TtL ITM Il9j| . 
The corresponding scaling laws are known from litera- 
ture [23, 123, 13 ■ The scaling equation for ground-state 
energy per site e = Eq /N of the bcc lattice is 

e(L) = e{oo) + A^L'^ + ... (7) 

and for the order parameter 

m'^ {L)^m^ {00) + B2L-^ + ... (8) 

with L = 



w -0.9 
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LSWT bcc 
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extrapolation bcc 1- 
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FIG. 6: LSWT data for the ground-state energies -Bo of the 
infinite bcc lattice (dashed line, LSWT bcc) and square lattice 
(dotted line, LSWT sq) as well as the extrapolated ED data 
for _Eo (solid line with data points) versus p. 
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with Wk = 

m = {Sf) is 



The sublattice magnetization where 



m 



N ^ 



1 

2 2tJk 



(15) 



a. LSWT for large .72 

The classical ground state of the Ji — J2 model on the 
bcc lattice for large J2 consists of two interpenetrating 
Neel states each living on the initial sublattices A and 
B. The two Neel states are energetically decoupled, i.e. 
the angle 9 between the staggered magnetization on A 
and B is arbitrary for classical spins. For the quantum 
model we start with arbitrary 9 and use as quantiza- 
tion axis the local orientation of the spins in the classical 
ground state. The further procedure is the same as for 
small J2 but the bosonic Hamiltonian now contains the 
angle 9. By means of the Hellmann-Feynman theorem 
[^] {dH{9)/d9) = dE{9)/de it can be easily found that 
in the quantum model the coUinear state (6' = or tt) 
has lowest energy. This lifting of the continuous degener- 
acy of the classical ground state by quantum fluctuations 

order from disorder effect ') is also found for the square 
lattice j26j|. For 9 — Q the bosonic Hamiltonian reads 



H = -6A^ps2 + ^ j A(at flk 6|,6k) (16) 
+ [Ck(6kak ^ ^k«^k) ^ Sk(aka_k + ^k^-k) + /i.c.jj 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



LSWT bcc 

extrapolation bcc — ' — 
LSWTsq 









A — 6sp 
Sk = 3sp72k 

Ck = 4s(cosA:j. cos/cy COS A:^ 
-|- i sin sin ky sin k^)- 



Then the ground-state energy per site is 
1 



with the modes 



N ^ 

k 



V 2 



^2k 

2 



A 



(17) 
(18) 
(19) 



(20) 



t^ik 



= ylA^-ABl 



-Fk ; W2k = _ 4^2 _ (21) 



and the function 



C*^ f - 8ABi,{Cl + 



(22) 



+ iCi,C^{A' + 4Bl)y^\ 
The sublattice magnetization is written as 



-E- 



^(k,c^ik) 



2tJlk(w2j^ - L 

N 2^ 



(23) 



N ^ 2tj2k(t^ 



2 

2k 



with 



0.2 



0.4 



0.6 



0.8 



D(k, c^k) - -ojI + AujI - (4^2 - A^)uo^ (24) 
+^(4^2 + 2CkQ) - 2Bk(C^ + C^) - A^. 

The results of LSWT and the finite-size extrapolation 
are shown in Figs.|Sl and [3 For the limits J2 = and 
Ji = our LSWT results are in agreement with data 
for the bcc and the sc lattice given in [23|. Both meth- 
ods yield similar results. For the ground-state energy we 
have a good quantitative agreement. Being in the size of 
the data points, the standard deviation of extrapolated 
ground-state energy is not shown. For the order parame- 
ter the finite-size effects are stronger and the agreement is 
only qualitative. Both methods suggest a first-order tran- 
sition for the spin-i J1-J2 model on the bcc lattice. The 
transition point obtained from the ED data is J2 ~ 0.7Ji 
(i.e. p « 0.52) while the LSWT becomes instable at the 
classical transition point. 



IV. CONCLUSION 



FIG. 7: LSWT data for the sublattice magnetizations m of the 
infinite bcc lattice (dashed line, LSWT bcc) and square lattice 
(dotted line, LSWT sq) as well as extrapolated ED data for 
^3m^(cxD) and ^3771^(00) (solid line with data points) vs. p. 
The error bars indicate the standard deviation of the finite- 
size extrapolation. 



We have presented spin-wave and exact diagonaliza- 
tion results for the spin-i Ji- J2 model on the bcc lattice 
and compare them with those for the square lattice. In 
general, we observe that the physics for the 3d quantum 
model is closer to classical behavior since quantum fluc- 
tuations and finite-size corrections become less important 
for higher coordination number and larger dimension. 
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We are not sure whether the increase of the magne- 
tization ttIq, approaching the transition point from the 
right (shown in Fig.|Z|) is a real effect. A possible physical 
origin for an increase may be a stronger coupling of the 
Neel ordered subsystems A and B due to larger quantum 
fluctuations or finite-size effects that become more 
important in the region of strong frustration. From the 
data for the ground-state energy, the Marshall-Peierls 
sign rule, the sublattice magnetizations and the spin 
gap we conclude that the increase from dimension d=2 
to d=3 changes the physical properties basically. The 
good agreement with the spin wave results support 
this conclusion. Contrary to the 2d model, where the 
quantum fluctuations and frustration lead to a second 



order transition from the two-sublattice Neel state to 
a disordered spin-liquid like phase, in the 3d model 
we find no indications for a disordered ground-state 
phase. The quantum Ji — J2 model on the bcc lattice 
shows one transition of first order induced by strong 
frustration from the two-sublattice Neel state directly to 
the collinear state, where the transition takes place at 
ac = {J2/Ji)c « 0.7. 
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